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Abstract 

This talk is based on my work in collaboration with Thibault 
Damour since 1991. Unified theories, like superstrings, predict the 
existence of scalar partners to the graviton. Such theories of grav- 
ity can be very close to general relativity in weak-field conditions 
(solar-system experiments), but can deviate significantly from it in 
the strong-field regime (near compact bodies, like neutron stars). Bi- 
nary pulsars are thus the best tools available for testing these theories. 
This talk presents the four main binary-pulsar experiments, and dis- 
cusses the constraints they impose on a generic class of tensor-scalar 
theories. It is shown notably that they rule out some models which are 
strictly indistinguishable from general relativity in the solar system. 
This illustrates the qualitative difference between binary-pulsar and 
solar-system tests of relativistic gravity. 
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1 Introduction 



The usual meaning of "testing a theory" is rather negative: One compares 
its predictions with experimental data, and a single inconsistency suffices to 
rule it out. On the other hand, it is difficult to determine what features of 
the theory are correct when it passes a given test. In order to extract some 
positive information from experiment, it is useful to embed the theory into a 
class of alternatives. Indeed, by contrasting their predictions, it is easier to 
understand in what way they differ, and to determine the common features 
which make them pass or not the available tests. Moreover, this approach 
can suggest new experiments to test the other features of the theories. 

The best known example of such an embedding of general relativity into 
a space of alternatives is the so-called Parametrized Post-Newtonian (PPN) 
formalism, which is extremely useful for studying gravity in weak-field condi- 
tions, at order 1/c 2 with respect to the Newtonian interaction. The original 
idea was formulated by Eddington who wrote the usual Schwarzschild 
metric in isotropic coordinates, but introduced some phenomenological pa- 
rameters /3 PPN , 7 PPN , in front of the different powers of the dimensionless 
ratio Gm/rc 2 : 
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General relativity, which corresponds to /5 PPN = 7 PPN = 1, is thus embed- 
ded into a two-dimensional space of theories parametrized by all real values 
of /? PPN , 7 PPN . [The third parameter that one may introduce in front of 
2Gm/rc 2 in goo can be reabsorbed in the definition of the mass m.] 

The constraints imposed in this space by solar-system experiments are 
displayed in Fig. |l], and give the following la limits on the Eddington pa- 
rameters: 

|/3 PPN -1| < 6xl(T 4 , (2a) 
| 7 PPN -1| < 2xl(T 3 . (2b) 

The PPN formalism has been further developed by Schiff, Baierlin, Nordtvedt 
and Will to describe any possible relativistic theory of gravity at order 1/c 2 . 
In particular, Will and Nordtvedt 0] introduced up to 8 extra parameters 
(besides /3 PPN and 7 PPN ), each of them describing a particular violation of 
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Figure 1: Solar-system constraints on the PPN parameters. The widths of 
the strips have been enlarged by a factor 100. The allowed region is shaded. 

the symmetries of general relativity, like local Lorentz invariance, or the con- 
servation of energy and momentum. Since these 8 parameters do not have 
any really natural field-theoretic motivation (as opposed to /9 PPN and 7 PPN ; 
see below), we will not consider them any longer in this paper. Let us just 
mention that they are even more constrained than (/? PPN — 1) and (7 PPN — 1) 
by solar-system experiments. [See the contribution of J. Bell to the present 
Proceedings for a discussion of the tight bounds on some of these parame- 
ters imposed by binary-pulsar data.] In the 10-dimensional space of all these 
PPN parameters (/? PPN , 7 PPN , and the 8 others), only a tiny neighborhood of 
Einstein's theory is thus allowed by solar-system experiments: the intersec- 
tion of the three strips of Fig. [I] and a very thin 8-dimensional slice parallel 
to the plane of this Figure. One can therefore conclude that general relativ- 
ity is essentially the only theory which passes all these tests, and one may 
naturally ask the question: Is it worth testing it any further ? 

The reason why solar-system tests do not suffice is the extreme weakness 
of the gravitational field in these conditions. Indeed, the largest deviation 
from the flat metric is found at the surface of the Sun, and is proportional 
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to its gravitational binding energy [Gmj Rc 2 ) Q «2x 10~ 6 (where R denotes 
the radius of the considered body). In the vicinity of the Earth, the grav- 
itational field is of order (Gm/Rc 2 )® ~ 7 x 10 -10 . This explains why only 
the first terms of the expansion (JJ) are tested by solar-system experiments. 
Two theories which are extremely close in weak-field conditions can differ 
significantly in the strong-field regime. For instance, the typical self-energy 
of a neutron star is Gm/Rc 2 ~ 0.2, and therefore one cannot justify any more 
the PPN truncation of the theory at order 1/c 2 . [A more rigorous definition 
of the binding energy is — <91nm/<91nG. This expression takes its maximum 
value, 0.5, for black holes. The value ~ 0.2 found for neutron stars should 
therefore be understood as a rather large number.] Binary pulsars are thus 
ideal tools for testing relativistic theories of gravity in strong-field conditions. 

Before embedding general relativity into a class of contrasting alterna- 
tives, and comparing their predictions with experimental data, let us first 
describe the four main binary-pulsar tests presently available. 

2 Binary-pulsar tests 

The aim of this talk is not to explain what is a pulsar to specialists of the 
question. For our purpose, it is sufficient to note that an isolated pulsar is 
essentially a (very stable) clock. A binary pulsar (a pulsar and a companion 
orbiting around each other) is thus a moving clock, the best tool that one 
could dream of to test a relativistic theory. Indeed, the frequency of the 
pulses is modified by the motion of the pulsar (Doppler effect), and one can 
extract from the Table Of Arrivals many information concerning the orbit. 
For instance, the orbital period P& can be obtained from the time between two 
maxima of the pulse frequency. One can also measure several other Keplerian 
parameters, like the eccentricity e of the orbit, the angular position uo of the 
periastron, etc. 

In the case of PSR B 1913+16, which has been continuously observed 
since its discovery in 1974 ||, the data are so precise that one can even mea- 
sure three relativistic effects with great accuracy, (i) The redshift due to the 
companion)] oc Gtub/t abc 2 and the second-order Doppler effect oc v\/2c 2 
are combined in an observable which has been denoted 7Timmg- [The index 
"Timing" is written to avoid a confusion with the Eddington parameter 7 PPN 
introduced in Eq. (|T5|).] Since the Keplerian parameters and ui have been 
measured accurately during two decades, their time derivatives are also avail- 
able: (ii) uj gives the periastron advance (a relativistic effect of order v 2 /c 2 ), 

1 A denotes the pulsar, B the companion, and tab the distance between them. 
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Figure 2: General relativity passes the (7Timing-^--ffc)i9i3+i6 test. 

and (iii) the variation of the orbital period, Pb, can be interpreted as a conse- 
quence of the energy loss due to the emission of gravitational waves (an effect 
of order v 5 /c 5 in general relativity, but generically of order v 3 /c 3 in alterna- 
tive theories; see below). The three "post-Keplerian" observables 7Timin g , 
tl), Pb can thus be compared with the predictions of a given theory, which 
depend on the unknown masses m A , m B of the pulsar and its companion. 
However, 3 observables minus 2 unknown quantities is still 1 test. The equa- 
tions lT iming {m A ,m B ) = 7ThLin g > u th {m A , m B ) = ^ obs , P b th {m A ,m B ) = P b obs , 
define three curves (in fact three strips) in the two-dimensional plane of the 
masses (m^We). If the three strips meet in a small region, there exists a 
pair of masses (m A ,m B ) which is consistent with all three observables, and 
therefore the theory is consistent with the binary-pulsar data. If they do not 
meet, the theory is ruled out. Figure || displays these strips in the case of 
general relativity, which passes the test with flying colors. [We will see below 
that some other theories can also pass this test.] 

The binary pulsar PSR B 1534+12 has been observed only since 1991 
[|j] but it is much closer to the Earth than PSR B 1913+16, and three 
post-Keplerian observables have already been measured with good precision: 
7Timing, and a new parameter denoted s. It is involved in the shape of the 
Shapiro time delay (an effect oc 1/c 3 due to the propagation of light in the 
curved spacetime around the companion), and it can be interpreted as the 
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Figure 3: General relativity passes the (jTimmg-u-s) 1534+12 test. 

sine s = sinz of the angle between the orbit and the plane of the sky. [The 
range r of this Shapiro time delay is also measured but with less precision.] 
Here again, the three strips "predictions* 11 (m^, tub) = observed values" can 
be plotted for a given theory, and if they meet each other, the test is passed. 
Figure ^ displays the case of general relativity, which passes the test at the 
la level §. 

As shown in Sec. 3 below, generic theories of gravity predict a large 
dipolar emission of gravitational waves (of order v 3 /c 3 ) when the masses of 
the pulsar and its companion are very different, whereas the prediction of 
general relativity starts at the much weaker quadrupolar order (oc t> 5 /c 5 ). 
Several dissymmetrical systems, like the neutron star-white dwarf binary 
PSR B 0655+64, happen to have very small observed values of Pb, consistent 
with general relativity but not with a typical dipolar radiation. This is the 
third binary-pulsar test that we will use to constrain the space of gravity 
theories. 

We will also see in Sec. 3 below that in generic theories of gravity, the 
acceleration of a neutron star towards the center of the Galaxy is not the same 
as the acceleration of a white dwarf. This violation of the strong equivalence 
principle causes a "gravitational Stark effect" on the orbit of a neutron star- 
white dwarf system: Its periastron is polarized towards the center of the 
Galaxy. [This is similar to the effects discussed in J. Bell's contribution to 
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the present Proceedings.] More precisely, the eccentricity vector e of the 
orbit is the sum of a fixed vector directed towards the Galaxy center 
(proportional to the difference of the accelerations of the bodies), and of 
a rotating vector e R (t) corresponding to the usual periastron advance at 
angular velocity Cjr. Several dissymmetrical systems of this kind (such as 
PSRs 1855+09, 1953+29, 1800—27) happen to have a very small eccentricity. 
The only explanation would be that the rotating vector e^(t) is precisely 
canceling the fixed contribution at the time of our observation: ep + 
en(t) ~ 0. However, this is very improbable, and one can use a statistical 
argument to constrain the space of theories ||. Moreover, by considering 
several such systems, the probability that they have simultaneously a small 
eccentricity is the product of the already small individual probabilities. This 
idea has been used in f7j to derive a very tight bound on the difference of 
the accelerations of the bodies. This is the fourth binary-pulsar test that 
we will use in the following. Of course, general relativity passes this test, 
since it does satisfy the strong equivalence principle (universality of free fall 
of self-gravitating objects). 

These four tests are presently the most precise of all those which are 
available. It should be noted that many other tests are a priori possible: 
Damour and Taylor || have shown that 15 tests are in principle possible for 
each binary pulsar, if the pulses are measured precisely enough. 

3 Tensor-scalar theories of gravity 
3.1 Introduction and action 

We saw in the previous section that several tests of gravity can be performed 
in the strong-field regime, and that general relativity passes all of them. As 
discussed in Sec. 1, our aim is now to embed Einstein's theory into a class of 
alternatives, in order to determine what features have been tested, and what 
can be further tested. A generalization of the PPN formalism to all orders 
in l/c n would need an infinite number of parameters [cf. Eq. (ID)]. We will 
instead focus on the most natural class of alternatives to general relativity: 
"tensor-scalar" theories, in which gravity is mediated by a tensor field {g^v) 
and one or several scalar fields ((p) . Here are the main reasons why this class 
is privileged, (i) Scalar partners to the graviton arise naturally in theoreti- 
cal attempts at quantizing gravity or at unifying it with other interactions 
(superstrings, Kaluza-Klein). (ii) They are the only consistent massless field 
theories able to satisfy the weak equivalence principle (universality of free fall 
of laboratory-size objects), (iii) They are the only known theories satisfying 
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"extended Lorentz invariance" 0, i.e., such that the physics of subsystems, 
influenced by external masses, exhibit Lorentz invariance. (iv) They explain 
the key role played by /3 PPN and 7 PPN in the PPN formalism (the extra 
8 parameters quoted in the Introduction vanish identically in tensor-scalar 
theories), (v) They are general enough to describe many different deviations 
from general relativity, but simple enough for their predictions to be fully 
worked out (T^ . 



Like in general relativity, the action of matter is given by a functional 
S m [ipm,g^i,} of some matter fields ip m (including gauge bosons) and one 
second-rank symmetric tensorQ g^ v . The difference with general relativity 
lies in the kinetic term of g^ u . Instead of being a pure spin-2 field, it is here 
a mixing of spin-2 and spin-0 excitations. More precisely, it can be written 
as 9fiv — exp[2a(</?)](7 M1/ , where a(y?) is a function of a scalar field <p, and g^ u 
is the Einstein (spin 2) metric. The action of the theory reads thus 



3 

S = I d A x^ (R - 2g» v d^d v ip) + S m ty m , e 2a ^g, v ] . (3) 

[Our signature is -+++, R is the scalar curvature of g^ v , and g its determi- 
nant.] 

Our discussion will now be focused on the function a((p), which charac- 
terizes the coupling of matter to the scalar field. It will be convenient to 
expand it around the background value <£> of the scalar field (i.e., its value 
far from any massive body): 

a((p) = a (y? - fo) + ^o(f ~ <Po) 2 + ~ fo) 3 + ••• , (4) 

where a , /3 , /3' , ... are constants defining the theory. General relativity 
corresponds to a vanishing function a(<£>) = 0, and Jordan-Fierz-Brans-Dicke 
theory to a linear function a(y?) = a (<^ — y^o)- We will see in Sec. 5 below 
that interesting strong-field effects occur when /3 ^ 0, i.e., when a(y?) has a 
nonvanishing curvature. 

3.2 Weak-field constraints 

Before studying the behavior of these theories in strong-field conditions, it 
is necessary to take into account the solar-system constraints (0). A simple 
diagrammatic argument |TT| allows us to derive the expressions of the effec- 



tive gravitational constant between two bodies, and of the Eddington PPN 



2 To simplify, we will consider here only theories which satisfy exactly the weak equiva- 
lence principle, and we will restrict our discussion to a single scalar field except in Sec. 4. 
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parameters in tensor-scalar theories: 



G 



cff 



7 



PPN 



PPN 



G(l+a 2 ) , 
-2ag/(l + a 
l_aoA)Oo 
2(l + ag 



2\2 



(5a) 
(5b) 

(5c) 



[The factor «q comes from the exchange of a scalar particle between two 
bodies, whereas a j3 a comes from a scalar exchange between three bodies.] 
The bounds (§) can therefore be rewritten as 



\alP 



< 
< 



1-3 



10" 

1.2 x 10~ 3 



(6a) 
(6b) 



The first equation tells us that the slope of the function a(ip) cannot be 
too large: The scalar field is linearly weakly coupled to matter. The second 
equation does not tell us much, since we already know that is small. In 
particular, it does not tell us if (3q is positive [a((p) convex] or negative [a(ip) 
concave] . 

The same diagrammatic argument can also be used to show that any 
deviation from general relativity at order 1 / c" (n > 2) involves at least two 
factors «o> an d has the schematic form 



deviation from G.R. = a x 



, Gm , fGm\ 2 



Rc 2 



\Rc 2 



(7) 



where Gm/Rc 2 is the compactness of a body, and Ao, Ai, ... are constants 
built from the coefficients ao, fio, /3' , . . . of expansion (§]). Since «q is known 
to be small, we thus expect the theory to be close to general relativity at 
any order. [We do not wish to consider models involving unnaturally large 
dimensionless numbers in the expansions (B) or (0).] However, in two differ- 
ent cases that will be discussed in Sections 4 and 5, the theory can exhibit 
significant strong-field deviations from general relativity: (i) If the theory 
involves more than one scalar field, Eq. (|6~aD does not necessarily imply that 
the slope of a(<p) is small, (ii) Some nonperturbative effects can develop in 
strong-field conditions, and the sum of the series in the square brackets of 
Eq. (0) can be large enough to compensate even a vanishingly small Oq. 



3.3 Strong-field predictions 

The predictions of tensor-scalar theories in strong-field conditions have been 
derived in ]T0| . They mimic the weak-field predictions with the important 
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difference that the constants aco, Po must be replaced by body-dependent 
parameters a a = dlnrriA/dipo, Pa = <9au/cVo (and similarly for the com- 
panion B) . These parameters can be interpreted essentially as the slope and 
the curvature of a((p) at the center of body A (or body B). [In the weak-field 
regime, one has ip ~ ip , therefore ~ a , Pa ~ Po-] In particular, the 
effective gravitational constant between two self-gravitating bodies A and B 
reads 

Gf B = G(l + a A a B ) , (8) 

instead of (f)a|). The acceleration of a neutron star A towards the center C of 
the Galaxy is thus proportional to (1 + aA®c), whereas a white dwarf B is 
accelerated with a force oc (1 + asctc)- Since 7^ «b in general, there is a 
violation of the strong equivalence principle which causes the "gravitational 
Stark effect" discussed in Sec. 2. 

The strong-field analogues of 7 PPN and P PPN are given by formulas similar 
to (|5T5|), (|5cl), but «q is replaced by oia&b and a p ao by a combination of 



o.aPbo<-a and ocbPa^b- The prediction for the periastron advance Co can thus 
be written straightforwardly. 

The expression of the observable parameter 7Timin g involves again the 
body-dependent parameters a a, o>b, but also a subtle contribution propor- 
tional to x d In Ia/O^q, where I a is the inertia moment of the pulsar. 
This term is due to the modification of the equilibrium configuration of the 
pulsar due to the presence of its companion at a varying distance. We have 



shown in [0 how to compute this effect, which happens to be particularly 
large in some models (see Sec. 5 below). 

The energy flux carried out by gravitational waves has been computed in 



It is of the form 



l Quadrupole 



Energy flux = < h O 

c 5 




helicity 2 

{Monopole Dipole Quadrupole / 1 \ 1 
~n + —^ + 75 +0 \7)( ■ ' ;<)) 



c c° c° \c 



helicity 

The first curly brackets contain the prediction of general relativity. The sec- 
ond ones contain the extra terms predicted in tensor-scalar theories. The 
powers of 1/c give the orders of magnitude of the different terms. In partic- 
ular, the monopolar and dipolar helicity-0 waves are generically expected to 
be much larger that the usual quadrupole of general relativity. However, the 
scalar monopole has the form 

Monopole _ G \ d(m A a A ) + d{m B a B ) + ^ / 1\ j 2 ^ 



c 



\ dt 
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and it reduces to order 0(l/c 5 ) if the stars A and B are at equilibrium 
[dtirriAOiA) = 0], which is the case for all binary pulsars quoted in Sec. 2. [It 
should be noted, however, that this monopole would be huge in the case of 
a collapsing star, for instance.] The dipole has the form 

and is usually much larger that a quadrupole of order 1/c 5 (see the third test 
discussed in Sec. 2). However, when the two stars A and B are very similar 
(e.g. two neutron stars), one has ~ and this dipolar contribution 
almost vanishes. [A dipole is a vector in space; two identical stars do not 
define a preferred orientation.] 



4 Tensor- mult i- scalar theories 

In order to satisfy the weak-field constraints (|) but still predict significant 
deviations from general relativity in the strong-field regime, the first possi- 
bility is to consider tensor-scalar theories involving at least two scalar fields 
[|l(J . Indeed, there can exist an exact compensation between the two fields in 
the solar system, although both of them can be strongly coupled to matter. 
If the kinetic terms of the scalar fields read — (c^i) 2 + (<9 M <^2) 2 , Eq. (|6a|) 
becomes \a\ — a\\ < 10~ 3 , and none of the coupling constants a\, a<i needs 
to be small. However, one of the fields (here (^2) must carry negative energy 
for this compensation to occur. Therefore, these tensor-bi-scalar theories 
can be considered only as phenomenological models, useful as contrasting 
alternatives to general relativity but with no fundamental significance. 

We have constructed in Ref. [|Hj] the simplest tensor-bi-scalar model which 
has the following properties: (i) It has the same post-Newtonian limit as 
general relativity (/5 PPN = 7 PPN = 1), and therefore passes all solar-system 
tests, (ii) It does not predict any dipolar radiation oc 1/c 3 [Wl, WB, (a a — 
ocb) 2 = 0], and therefore passes the third binary-pulsar test discussed in 
Sec. 2. Moreover, it depends on two parameters, /?', f3", and general relativity 
corresponds to (3' = (3" = 0. Figure || displays the constraints imposed by 
the three other binary-pulsar tests of Sec. 2 in the plane of the parameters 
(/?',/?"). The theories passing the 1913+16 test are inside the long strip 
plotted in solid lines. Note that theories which are very different from general 
relativity can pass this test. For instance, Fig. [5] displays the mass plane 
(myijine) for the (fine-tuned) model f3' = 8, (3" = 69. The three strips are 
significantly different from those of Fig. 0, but they still meet each other in 
a small region [corresponding to values of the masses m^, m# different from 
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Figure 4: Constraints imposed by the binary-pulsar tests of Sec. 2 on the 
tensor-bi-scalar model of Sec. 4. The two dotted strips illustrate how the 
precision of the "Stark" test is increased when several binary pulsars are 
considered simultaneously. The region allowed by all tests is the small shaded 
diamond around general relativity (/?' = f3" = 0). 

those found in general relativity]. To illustrate how much the theory differs 
from general relativity, let us just mention that the effective gravitational 
constant G e ^ B between the pulsar and its companion is 1.7 times larger than 
the bare Newtonian constant G. We have thus exhibited a model which 
deviates by 70 % from Einstein's theory, but passes (i) all solar-system tests, 
(ii) the "no-dipolar-radiation" test of PSR 0655+64, and (iii) the "7Timin g - 
iu-Pb" test of PSR 1913+16. Before our work, this 1913+16 test was usually 
considered as enough to rule out any theory but general relativity. We have 
proven that other binary-pulsar tests are also necessary. In particular, Fig. ^ 
shows clearly that the "7Timing-^-s" test of PSR 1534+12 and the "Stark" 
test complement it usefully. For instance, the model of Fig. |5] is easily ruled 
out by the 1534+12 test: the 7Timin g and s curves do not even meet each 
other (so that the observable u> is not even useful here). 
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Figure 5: The tensor-bi-scalar model (3' = 8, (3" = 69, passes the (7Timm g - 
Co-Pb) 1913+16 test, although the three curves are significantly different from 
those of Fig. ||. 

Thanks to the four binary-pulsar tests discussed in Sec. 2, this class of 
tensor-bi-scalar models is now essentially ruled out. We have achieved a 
similar result as in the weak-field regime of Fig. |l|: Only a tiny neighborhood 
of general relativity is still allowed. This is a much stronger result than just 
verifying that Einstein's theory passes these four tests. 

5 Nonperturbative strong-field effects 

We now discuss the second way to satisfy the constraints (^|) while predicting 
significant deviations from general relativity in the strong-field regime. As 
opposed to the models of the previous section, we consider here well-behaved 
theories, with only positive-energy excitations (of the type that is predicted 
by superstrings and extra-dimensional theories). To simplify, we will also 
restrict our discussion to the case of a single scalar field ip. 

The simplest tensor-scalar theory, Jordan-Fierz-Brans-Dicke theory, can- 
not give rise to nonperturbative strong-field effects for an obvious reason. It 
corresponds to a linear coupling function a(y?) = a ((p — <p ), and even if the 
field ip a at the center of body A is very different from the background <p , 
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one has anyway a a ~ a' ('Pa) = «o- Therefore, the deviations from general 
relativity, proportional to ~ a^, are constrained by the solar-system 

limit ( pip to be <^ 0.1 % even in the vicinity of neutron stars. 

On the contrary, if we consider a quadratic coupling function a(tp) = 
|/3o<£> 2 , the field equation for ip in a body of constant density p is of the form 
d 2 (rp)/dr 2 « /3 p • (^V 9 )- Therefore, the solution involves a sinh if /3 > 0, 
and a sin if /3 < 0. More precisely, one finds 



ola ~ <W cosh J 3 (3 Gm/ Re 2 if /3 > , (12a) 



a A w Qt / cos v /3|/3 |Gm/ J Rc 2 if /3 < . (12b) 

In the case of a convex coupling function a(tp) (i.e., P > 0), the deviations 
from general relativity are thus smaller in strong-field conditions that in the 
weak-field regime: a^ag < «q < 1CT 3 . On the other hand, a concave a(<p) 
can give rise to significant deviations: If /3 <; —4, the argument of the cosine 
function is close to tt/2 for a typical neutron star (Gm/Rc 2 ~ 0.2), and a A 
can thus be large even if a is vanishingly small. To understand intuitively 
what happens when a = strictly (i.e., when the theory is strictly equivalent 
to general relativity in weak-field conditions), it is instructive to compute the 
energy of a typical configuration of the scalar field, starting from a value <pa 
at the center of body A and tending towards as 1/r outside. One gets a 
result of the form 



Energy 



i(%>) 2 + pe^'A « mc 2 ( -Ml 



When (3q < 0, this is the sum of a parabola and a Gaussian, and if the com- 
pactness Gm/Rc 2 is large enough, the function Eneigj(pA) has the shape 
of a Mexican hat; the value <pa = now corresponds to a local maximum 
of the energy. It is therefore energetically favorable for the star to create 
a nonvanishing scalar field pa, and thereby a nonvanishing "scalar charge" 
ola ~ PoVa- This phenomenon is analogous to the spontaneous magnetiza- 
tion of ferromagnets. 

We have verified the above heuristic arguments by explicit numerical cal- 



culations [pL3| , taking into account the coupled differential equations of the 
metric and the scalar field, and using a realistic equation of state to de- 
scribe nuclear matter inside a neutron star. We found that there is indeed a 
"spontaneous scalarization" above a critical mass, whose value depends on 
Pq. Figure |] displays the scalar charge a a for the model Pq = —6. Note that 
the deviations from general relativity are of order olaolb ~ 35 % for a wide 
range of masses from w 1.25 m© to the maximum mass; therefore, no fine 
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Figure 6: Scalar charge a a versus baryonic mass tua, for the model a((p) = 
—3ip 2 (i.e., Pq = —6). The solid line corresponds to the maximum value 
of a allowed by solar-system experiments, and the dashed lines to a = 
("zero-mode"). The dotted lines correspond to unstable configurations of the 
star. 

tuning is necessary to get large deviations in a particular binary pulsar. Note 
also that the nonperturbative effects do not vanish with «o : Even if the the- 
ory is strictly equivalent to general relativity in the solar system, it deviates 
significantly from it near compact bodies. In fact, an even more surprising 
phenomenon occurs for the term a_e9 In 7^/9(^0 involved in the observable 
7Timin g (see Sec. 2): This term blows up as a — ► 0. In other words, a theory 
which is closer to general relativity in weak-field conditions predicts larger 
deviations in the strong-field regime ! 

The "7Timin g -^-A" test of PSR 1913+16 is displayed in Fig. for the 
model (3q = —6 and the maximum value of «o allowed by solar-system ex- 
periments. [The "7Timing-^-s" test of PSR 1534+12 gives curves similar to 
those of Fig. |7| for the first two observables, while the s strip is only slightly 
deviated from that of Fig. ||.] The great deformation of the Pb curve, as 
compared to the general relativistic prediction, Fig. |2], is due to the emission 
of dipolar waves in tensor-scalar theories. The fact that this dipolar radia- 
tion vanishes on the diagonal = mg explains the shape of this curve. As 
expected, the 7Timin g curve is also very deformed because of the contribution 
asd In I^/dtpo. When f3 is not too negative (e.g. @ ~ —4), a smaller value 
of a allows the test to be passed: the three curves finally meet in one point. 
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Figure 7: The model a(ip) = —3f 2 does not pass the (7Timing-^-A)i9i3+i6 
test. 



On the contrary, when A < —5, we find that the test is never passed even 
for a vanishingly small cto (because the term aedlni^/cfyjo blows up). In 
other words, this binary pulsar rules out all the theories A < —5, a = 0, 
although they are strictly equivalent to general relativity in weak-field con- 
ditions. This illustrates the qualitative difference between binary-pulsar and 
solar-system tests. 

Generic tensor-scalar theories can be parametrized by the first two deriva- 
tives, «o and Po, of their coupling function a((p), cf. Eq. (||). It is instructive 
to plot the constraints imposed by all kinds of tests in the plane (ao, AO- 
Figure H] shows that solar-system experiments do not constrain at all the 
curvature A) °f ^V 9 ) if its slope ao is small enough. On the contrary, binary 
pulsars impose A > — 5, independently of ato- Using Eqs. (|5b|) , (|5c|), this 
bound can be expressed in terms of the Eddington parameters: 



A 



PPN 



1 



7 



PPN 



< 1.3 . 



(14) 



The singular (0/0) nature of this ratio vividly expresses why such a conclusion 
could not be obtained in weak-field experiments. 



Recent cosmo logical studies, notably [14], have shown that theories with 
a positive A are easily consistent with observational data, whereas some 
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general 
relativity 

Figure 8: Constraints imposed by solar-system and binary-pulsar experi- 
ments in the plane (an, A))- m view of the reflection symmetry a — > — a , 
only the upper half plane is plotted. The allowed regions are below and on 
the right of the different curves. The shaded region is allowed by all the tests. 

fine-tuning would be required if /3n < 0. It is fortunate that binary pulsars 
precisely privilege the positive values of this parameter. 

6 Conclusions 

Tensor-scalar theories of gravity are the most natural alternatives to general 
relativity. They are useful as contrasting alternatives, and can suggest new 
experimental tests. For instance, the tensor-bi-scalar model of Sec. 4 proved 
that a single binary-pulsar test does not suffice. Well-behaved tensor-scalar 
theories (with no negative energy, no large dimensionless parameters, and no 
fine tuning) can develop nonperturbative strong-field effects analogous to the 
spontaneous magnetization of ferromagnets. Their study illustrates the quali- 
tative difference between binary-pulsar and solar-system experiments: binary 
pulsars have the capability of testing theories which are strictly equivalent 
to general relativity in the solar system. 
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